Coherent electronic transport in a multimode quantum channel 
with Gaussian-type scatterers 
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Coherent electron transport through a quantum channel in the presence of a general extended 
scattering potential is investigated using a T-matrix Lippmann-Schwinger approach. The formalism 
is applied to a quantum wire with Gaussian type scattering potentials, which can be used to model 
a single impurity, a quantum dot or more complicated structures in the wire. The well known dips 
in the conductance in the presence of attractive impurities is reproduced. A resonant transmission 
peak in the conductance is seen as the energy of the incident electron coincides with an energy level 
in the quantum dot. The conductance through a quantum wire in the presence of an asymmetric 
potential are also shown. In the case of a narrow potential parallel to the wire we find that two dips 
appear in the same subband which we ascribe to two quasi bound states originating from the next 
evanescent mode. 

PACS numbers: 72.10.Bg,73.63.Nm,73.63.Kv 



I. INTRODUCTION 

In studies of electronic transport in mesoscopic or 
nanoscale quasi-one-dimensional systems it has strikingly 
been found that the conductance manifests quantiza- 
tion^i^i^ when the electronic phase-coherent length 
is greater than the system size. The transport prop- 
erties of these systems have been successfully described 
within the Landauer-Biittiker frameworki^iSiiSiSiifl Later 
on, the influence of a single impurity on the conduction in 
quantum channels has attracted a great deal of attention 
since the impurity inside or near the conducting channel 
may destroy the conduction quantization, as has been 
demonstrated theoretically Iiii2ii3ii4,i5,i6,i7,i8,i9 
perimentallyj^U 

The impurities are usually assumed to be zero- 
range, i.e. of a delta-function type, in theoretical 
considerationsiii*2iiS^ However, a real physical impurity 
must be of finite range and therefore modeled by an ex- 
tended potential. In this paper we discuss coherent trans- 
port properties of quantum channels in the presence of 
extended scattering potentials. If the potentials are very 
extended they can describe the effects of a tunable central 
gateiS2i 

The method for the coherent electron transport 
through the quantum channel employs the Lippmann- 
Schwinger (LS) formalism. Its relation to the T-matrix 
and the Landauer formula in terms of elements of the 
scattering matrix is discussed in Sec. ^ The conduc- 
tance of a quantum channel in the presence of extended 
Gaussian type scattering potentials obtained by numeri- 
cal calculations is given in Sec. lIIII This section consists 
of three parts. First we discuss the scattering by a sin- 
gle Gaussian potential. The physics of the transport is 
examined more closely by visualizing the most impor- 
tant processes by help of the LS equation. In the second 
part we model a quantum dot embedded in the wire by 



a combination of two Gaussian potentials, and use that 
to obtain the conductance of such a setup. The effect 
of the shape of an extended scattering potential on the 
conductance is considered in the last part. In Sec. II VI we 
summarize and discuss the main results of the paper. 



II. SCATTERING BY A POTENTIAL 

We consider a quantum wire connected adiabatically 
to reservoirs as schematically shown in Fig. The 
uniformity of the wire is broken by a scattering po- 
tential of finite extent in the wire. The wave func- 
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Figure 1: Schematic view of the system. An incoming wave 
is partly transmitted and partly reflected by the flnite range 
scattering potential V. 



tion of a single-electron state with energy E, in the 
quasi-one-dimensional quantum wire is described by the 
Schrodinger equation 
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+ Vc{x) + V^(r) ) i;E{v) = E^e{v). 

(1) 

The electron is confined in the x-direction by the con- 
finement potential Vc{x) but is free to propagate in the 
z-direction. V{y) is a finite range scattering potential. 
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The transverse modes 



2m dx"^ 



+ Vcix) Xn{x) = SnXnix), 



(2) 



are assumed to be known. 

Outside the range of the scattering potential an elec- 
tron with energy is in a linear combination of the eigen- 
functions of Hq = H — V{y). These eigenfunctions are 
referred to as modes and can be written 
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Xn{x), 



(3) 



where E = 1i'k'^[E)/2m + e„ and the modes are nor- 
malized to carry unit probability current The -I- (-) 
refers to waves incident from the left (right). If e„ < E, 
kn{E) is real and (j)^^ are propagating waves (see Fig.|2l. 
If on the other hand e„ > E, kn{E) is purely imaginary 
and the eigenfunctions must be exponentially decaying, 
evanescent modes. 




Figure 2: (Color online). The energy levels of a quantum 
wire with a hard wall confinement. For a given energy E 
(horizontal line) only a finite number of propagating modes 
are active (red markers), but higher lying evanescent modes 
can serve as intermediate states. This specific case is for a 
wire of width L — 4ao. Energy is scaled in Rydbergs (Ryd) 
and lengths in Bohr radii (ao), cf. Sec. IIIII 

Scattering states are special solutions of Eq. (QJ , where 
e.g. an electron in propagating mode n impinges on the 
scattering potential from the left and is backscattered 
or transmitted into other modes. In the far zone, outside 
the range of the scattering potential, the scattering states 
have the boundary condition 



m,prop 



tmn<l>mEi^)-' 



OO. 



(4) 

t„in is the probability amplitude for a particle in mode 
n in the left lead to scatter out in mode m in the right 
lead, rmn gives the probability amplitude for reflection 
into mode to. By restricting all the modes to carry unit 
current the transmission amplitudes coincide with the 



elements of the scattering matrixiiS Hence, the zero tem- 
perature linear response conductance of the system can 
be obtained by the Landauer formulaifl*2Si 



h 



(5) 



where the matrix of the transmission amplitudes t is eval- 
uated at the Fermi energy. 



A. Scattering States and the Lippmann-Schwinger 
Equation 

We expand the scattering states in the transverse 
modes 



i^nEi'^) = <^ms(^)Xm(a;), 



(6) 



where the n denotes the number of the incident mode (cf. 
Eq. Q). Introducing this expansion in the Schrodinger 
equation multiplying with Xm'i^)^ integrating over 
X and letting m ^ m' , one obtains a coupled mode 
equationiiiSS 

(7) 
(8) 

(9) 



where 

Vrr,m'{z) ^ J dx Xrn{x)V {r)Xm' (x) ■ 

Defining a mode Green's function as 

(^■^+kfAE))GnEiz,z')=Siz-z'), 



the solution to Eq. can be written in the form of an 
effective ID LS equation 



^rnEiz) = ^7nE{z) + ^ 

n 



E / dz'g^Eiz^z') 

m' 

X V^^,{z')ip'^,E{z'), 



(10) 



where (f'^Eiz) = ^nm exp{ik,n{E)z)/ ^k„i{E). Note that 
the solutions of the LS equation ifTfljl obey the same nor- 
maHzation as the incident mode. Inserting the explicit 
form of the Green's function. 



G°nE{z,z') 



^ik„{E)\z-z'\ 



2kniE) 



(11) 



in Eq. IjlOII and taking the limit z ^ oo, one obtains by 
comparison with Eq. 1^1 the probability amplitudes for 
forward scattering 



-Y^jdz'^ 



-ik^(E)z' 



xV.,m'{z')^l,.{z'). 



(12) 
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B. Scattering States and the T-matrix 

Having obtained the transmission amplitudes in terms 
of the configuration space wave function we now find 
them in terms of matrix elements of a transition opera- 
tor. Inserting the definition ||HJ of Vmm' into relation lfT2|l 
and using the expansion ijHJ one obtains 



(mfc™|]/|^+^), 



(13) 



where in general {r\nq) — {exp{iqz) / ^/\q\)xm{x) . Simi- 
larly, using the LS equation in the same expansion, 
we see that 



(14) 



where 



2 m 



Goir, r'; i?) = tt E xU^)G"raE{^, z')xm{x). (15) 



This is just the conventional LS equation and therefore 
by defining the T-matrix by the relation 



{mk'\f\nK) = (mfe'lFlV+s), 
it satisfies the operator LS equation 

f = v + vgQ[E)f. 



(16) 



(17) 



In the eigenfunction basis of i?o the LS equation ltT6| is 
transformed into2^ 



Tmn{k, kji) — Vmn{k, k^) H — o / I dq \q\ 



Vmiik,q)Ti„{q, k„) 
kf — q^ + iri 



(18) 



where the notation Vmi{k,q) = (mfc| V \lq) has been in- 
troduced. We have used that in our current normaliza- 
tion the unity operator is 

i = E (19) 



and that the Green's function is 

l-K SwSiq-q') 



{lq\go{E)\l'q') 



- ei + irj 



(20) 



The numerical solution of Eq. ifTsll is briefiy discussed in 
App.EI 

In the above notation the transmission amplitudes can 
be written 



tr. 



(21) 



Note that through Eq. Ijl8ll Tmn{km,kn) depends on 
T-matrix elements both on and off the energy shell 
(h'^k'^ /2m + En is equal to E on the energy shell but 
not off it) . 

We have thus managed to link the elements of the 
transmission matrix to matrix elements of the transition 
operator T. The result is intuitively appealing, transmis- 
sion from one channel to another is related to the transi- 
tion from an eigenstate in one channel to an eigenstate in 
another. The formalism is quite general with respect to 
the type of wire confinement and scattering potential of 
finite range. A change of either one requires only a new 
evaluation of matrix elements of the scattering potential 
and the energy spectrum of Hq. 



III. RESULTS 

In our model we use two potentials, one describing the 
confinement and the other representing the scattering po- 
tential. We use two different models for either one of 
these. The confining potential, characterizing the extent 
of the quantum channel, can be described by either a 
hard wall potential 



V'c,hard(a;) 



Q<x<L 
oo elsewhere. 



(22) 



with eigenvalues e„ = ?i^7r^(n/T)^/2m, rt = 1, 2, . . ., or 
parabolic walls 



K,paraboiic(a;) = \^mu?x^r. 



(23) 



with eigenvalues £„ = (n + l/2)?iti>, n = 0, 1, . . .. Note 
that the first mode in the hard wall confinement has n = 
1 but n = in the parabolic case. 

By a proper choice of model for the scattering potential 
one can study a multitude of different systems. In the fol- 
lowing subsections we will by use of Gaussian functions 
model a single impurity by a single Gaussian, a quantum 
dot embedded in the wire by a sum of two Gaussians of 
different size and finally study shape effects of the scat- 
tering potential by asymmetric Gaussians. Even though 
we have chosen to use only Gaussian potentials we stress 
that the formaHsm is general and can be used with more 
complicated potentials. By cleverly choosing the func- 
tional form of the scattering potential it can be adapted 
to describe many experimental setups. 

In this paper we scale all energies in Rydbergs (Ryd) 
and lengths in units of Bohr radii (ao). The calculations 
are independent of exact material constants but it can be 
useful to keep in mind that in GaAs oq — 9.79 nm and 
Ryd = 5.93 meV. 
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A. Single Gaussian Potential 

A single Gaussian potential 

V{r) = T/oe-"((="-^')'+"'). (24) 

can be used to model an impurity in the wire or the 
effect of a tunable central gatCf^ In this subsection and 
the ones that follow, the center (xi,0) of the potential 
is chosen to be in the middle of the wire (i.e. Xi — L/2 
and Xi = in the hard wall and parabolically confined 
quantum wires respectively) unless otherwise noted. 




E/Ryd 



Figure 3: (Color online). Conductance of a hard wall quan- 
tum wire in units of Go = 2e^//i as a function of energy in the 
presence of a single repulsive (dashed, Vo — 3 Ryd, a — 2ag^) 
and attractive (solid, Vo = —4 Ryd, a — ^a^^) Gaussian scat- 
tering potential. The inset shows the scattering potentials in 
a cross section along the middle of the wire (a; = L/2). The 
total number of modes is A'' = 8. The width of the wire is 
L = 4ao. 

The conductance of a hard wall quantum wire in the 
presence of such a Gaussian scatterer is shown in Fig. El 
In the case of a repulsive potential, the well known steps 
in the conductance are little smeared out, similar to what 
is seen in experiments^ One should remember though 
that our calculations are at zero temperature, so there 
is no temperature smearing. In experiments, the finite 
temperature can be a significant source of smearing of 
the steps. Additionally, in the case of an attractive po- 
tential, there are dips in the conductance right before 
the onset of the next mode. The dips can be under- 
stood from a simple coupled mode model to be due to 
backscattering by a quasi-bound state originating from 
an evanescent mode in the next energy subband i^^i^^ A 
more general explanation in terms of the symmetries of 
the scattering matrix also exist iSi This effect is a multi- 
mode effect that disappears if evanescent modes are not 
included in the calculations. Interestingly, there is no dip 
in the first subband. This can also be understood from 
the coupled mode model, in which coupling between the 



first propagating mode and the quasi-bound state in the 
first evanescent mode is given by Vi2{z). Due to symme- 
try this matrix element is zero, the potential being even 
and the wave functions for the first and second transverse 
modes being even and odd respectively, hence no dip. 
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Figure 4: (Color online). Conductance of a parabolic quan- 
tum wire in units of Go = 2e^//i as a function of energy 
in the presence of a single repulsive (dashed, Vo — 1.35 
Ryd, a = 0.96a^^) and attractive (solid, Vo ~ —1.35 Ryd, 
a = 0.96aQ ^) Gaussian scattering potential. The inset shows 
the scattering potentials in a cross section along the middle 
of the wire {x = 0) and the parabolic confinement potential 
{hiu = 1.01 Ryd). The total number of modes is A'^ = 9. 
The arrows point at the the energies at which the probability 
density is calculated in Figs. Eland El 

The conductance of a parabolic wall quantum wire in 
the presence of a single Gaussian scattering potential is 
shown in Fig.^ Due to different dependence of the trans- 
verse energy levels on the subband number (linear and 
square) the length of the plateaus increases in the case 
of hard wall potential while being of constant length for 
the paraboHc wires. Besides this, there is no qualitative 
difference between the two confinement potentials. 

The conductance of the hard wall quantum wires 
has been obtained by solving for the T-matrix through 
Eq. Ijl8|l . The conductance of the parabolic wall quantum 
wire has both been obtained by the T-matrix approach 
and by solving the real space LS Eq. itTTHl . In the latter 
case the real space wave function is a part of the solu- 
tion but has to be calculated separately in the T-matrix 
case (cf. App.EJ- Visualizing the wavefunctions can aid 
in discussing the scattering processes. Therefore, we now 
visually examine the system of the paraboHcally confined 
wire at the energies marked by an arrow in Fig. 31 making 
the physics of the process more clear. 

In Fig. we have plotted the probability density 
[■(/'^^(r)!^ of the scattering state for incoming waves with 
n = and n = 1, at the energy indicated by the left 
arrow in Fig. 31 

The incoming wave is partly refiected and partly trans- 
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Figure 5: (Color online). The probability density of the scat- 
tering states tp^^ in the parabolic quantum wire in the pres- 
ence of the single attractive Gaussian scattering potential of 
Fig. 2] The total energy of the incident particle is i5 = 2.02 
Ryd (cf. the left arrow in Fig. The incoming wave is in 
mode n = (above) and n = 1 (below). 



Figure 6: (Color online). The probability density of the scat- 
tering states V'^E in the parabolic quantum wire in the pres- 
ence of the single attractive Gaussian scattering potential of 
Fig. 2] The total energy of the incident particle is i5 = 2.50 
Ryd, coinciding with a dip in the conductance (cf. the right 
arrow in Fig.l^. The incoming wave is in mode n = (above) 
and n — 1 (below). 



mitted. To the left of the scattering center the incoming 
part and the reflected part interfere to form the beating 
pattern seen. To the right, only the transmitted part ex- 
ists and the probability density is independent of z. It is 
interesting to note that in general the transmitted part 
in the second subband is a linear combination of the first 
two modes. The two modes should then interfere to make 
a beating pattern also on the right hand side. However, 
as already mentioned, the coupHng between the first two 
modes is zero due to symmetry. Therefore an incident 
electron in the first mode cannot be transmitted into the 
second mode. Examining a similar picture in the third 
subband reveals interference between mode n — and 
mode n = 2 in the right hand side, in agreement with 
the above. The effect of the normal modes can also be 
seen in the nodal structure in the transverse direction. 
The second mode, n = 1, has a node in the middle of the 
wire while the first mode, n = 0, is nodeless. 

More interesting things are happening in the dip in 
the second subband. The upper panel of Fig. |H1 shows 
the probability density of the scattering state obtained 
when the incident wave is in mode n = 0. This mode is 
resonantly backscattered due to coupling to an evanes- 



cent state in mode n — 2. Mode n = 1 is unaffected be- 
cause of symmetry blocking (lower panel of Fig. . The 
evanescent state with the corresponding nodal structure 
is clearly seen in the total wavefunction as a locaHzed 
state around the scattering potential. A beating pattern 
to the left of the scattering center, similar to the one seen 
in Fig. El due to the interference of the incoming and re- 
flected wave is faintly seen. On the other side of the po- 
tential, however, no outgoing wave is seen in agreement 
with resonant backscattering. An incident wave in mode 
n = 1 can due to symmetry not couple to the evanes- 
cent state. A scattering state with such an incident wave 
is thus expected to be qualitatively the same as for an 
energy away from the dip in the same subband, as can 
be seen in the lower panel of Fig. The main differ- 
ence between the scattering states with incident wave in 
mode n = 1 in Fig. and Fig. is that the latter is at 
higher energy. The wave vector is thus higher, resulting 
in higher frequency in the oscillations of the probability 
density. The same is true when comparing the different 
frequencies of the oscillations in the two cases of Fig. El 
In the n = smaller part of the energy is in the 
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transverse mode thus increasing the wave vector of the 
wave in the propagating direction. 

Both of these effects, the steps and the dip, are well 
known in the literature. However, in most calculations 
the scattering potential is of zero-range delta function 
type i^^i^^'^^ Extended potentials are rare, an example 
being the rectangle potential of Bagwell.-^ An interest- 
ing difference can be seen in the total number of modes 
needed in the calculations. In our case, 8 modes give re- 
sults that change insignificantly if more modes are added. 
In the delta function case, up to 100 modes are needed, a 
clear signature of the singular nature of the delta function 
potential. 

B. A Quantum Dot Embedded in the Wire 
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Figure 7: (Color online). Conductance of a hard wall quan- 
tum wire in units of Go ~ 2e^ /h as a function of energy in 
the presence of a double Gaussian scattering potential with 
a varying depth of the well. The inset shows the scattering 
potentials, whose parameters in the solid (dashed) case are 
Vi = 5.21 (6.16) Ryd, V2 = 7.49 (9.70) Ryd, Qi = 0.52 (0.60) 
Sq ^ and Q2 ~ 1.52 (1.60) Hq ^, in the cross section x — L/2. 
The total number of modes is A'^ = 8. The width of the wire 
is L = 4ao. 

The double Gaussian 

V{r) ^Vie-"'^'-' -¥26-"^"-^ a2>ai, (25) 

can be used to model a quantum dot embedded in the 
wire. The conductance of a hard wall quantum wire with 
such a quantum dot is given in Fig. [3 Compared to the 
single impurity, a new effect appears. An enhanced con- 
ductance in the first subband along with similar weaker 
structures in higher bands. This is the well known effect 
of resonant tunneling. A remnant of a bound state in the 
well induces a resonance in the transmission, hence the 
resonance in the conductance. Examining the coupled 
mode equation 0, one notices that, ignoring the cou- 
pling between the modes, the transmission is governed 



Figure 8: (Color online). The solid line represents the double 
Gaussian potential of Fig. along with its first two effective 
potentials. 



by the effective potentials Vnn{z). In Fig. we have 
plotted the double Gaussian potential along with its first 
two effective potentials. We know from resonant tunnel- 
ing in ID, that the higher the probability for the electron 
to tunnel out of the well, the broader the resonance in 
the transmission. The difference in the sharpness of the 
tunneling resonances in different subbands can therefore 
be attributed to the different shapes of the effective po- 
tential. 




Figure 9: (Color online). Conductance of a parabolic wall 
quantum wire in units of Go — 2e^ /h as a function of energy 
in the presence of a double Gaussian scattering potential. The 
inset shows the scattering potential, whose parameters are 
Vi = 2.02 Ryd, V2 = 3.71 Ryd, ai = 0.29 a^^ and 02 = 0.96 
^, in the cross section a; = 0. The parabolic confinement 
potential {hco = 1.01 Ryd) is also shown. The total number of 
modes is A'^ = 9. The arrow points at the value of the energy 
at which the probability density in Fig. ED is calculated. 
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The conductance through a quantum dot in a paraboHc 
wire m Fig. H has the same pronounced resonance in the 
first subband but the resonances in the higher subbands 
are even more suppressed than in the hard wall wire. 
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Figure 10: (Color online). The probability density of the 
scattering states ipts ™ parabolic quantum wire in the 
presence of the double Gaussian scattering potential of Fig. El 
The total energy of the incident particle is E — 0.64 Ryd, 
coinciding with a resonance in the conductance (cf. the arrow 
in Fig.EJ. The incoming wave is in mode n = 0. 

The resonant tunneling is through the remnants of 
a bound state in the well. In contrast to the reso- 
nant backscattering, this quasi-bound state belongs to 
the propagating subband, not the next evanescent mode. 
The resonance is thus not a multimode phenomenon as 
the dip certainly is. We examine the probability density 
in the most pronounced resonance of Fig. |2l that is the 
one marked by an arrow. There is only one propagating 
mode and thus the only scattering state is the one with 
an incident wave in mode n = 0; the probability density 
of which is shown if Fig. The quasi-bound state in 
the middle of the wire is clearly seen. The absence of 
beating on both sides of the scattering center indicates 
that there is no reflected wave, i.e. the incoming wave is 
perfectly transmitted. This supports very well the sug- 
gestion that the resonance is due to resonant tunneling 
through a remnant of a bound state in the well. 



C. Shape Effects of the Scattering Potential 
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Figure 11: (Color online). Conductance of a hard wall quan- 
tum wire in units of Go = 2e^//i, as a function of energy. A 
single contour of the scattering potential, an asymmetric sin- 
gle Gaussian, is shown in the inset. The contour is chosen to 
be at the value of Vbe~^ = 0.368Vb. In the upper panel the 
potentials are repulsive (Vb — 3 Ryd) but attractive (Vo ~ —3 
Ryd) in the lower panel. In both cases the circularly symmet- 
ric potential has a = 2 a^^, the transverse one has ~ 0.10 
Oq ^ and Qz = 6.73 a^^, and the longitudinal has ax = 9.10 
0.11 Uq^. The parameters are chosen such that 



ag and 



the volume between the potential and the zero energy plane 
are equal. The width of the wire is L = 6 ao and the total 
number of modes is A*" = 8. The symmetric and longitudinal 
potentials have been shifted by 1 Go and 2 Go respectively, 
for clarity. 



So far, the scattering potentials under examination 
have all been circularly symmetric. Elongated extended 
scattering potentials, such as the asymmetric Gaussian 

y(r) = V^oe~"-("-"*''-"^"', (26) 

offer more variety. The conductance characteristics for 
the cases of elongated longitudinal, elongated transverse, 
and circularly symmetric potential profiles are under in- 
vestigation, as shown in Fig. El The parameters of the 
potentials have been chosen such that their volume in 



the wire is the same in all cases. This is done so that the 
strength of the scatterers is comparable in all cases. We 
have used both attractive and repulsive potentials. 

We begin by discussing the conductance of the repul- 
sive potentials in the upper panel of Fig.^] The conduc- 
tance of the transverse barrier is lower than its circularly 
symmetric counterpart. The conductance of the longitu- 
dinal barrier has a more pronounced exponential growth 
and is thus lower than the conductance of the circularly 
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symmetric barrier at the low energy part of each subband 
while being higher in the high energy end. Both the sym- 
metric and the longitudinal barrier are mainly confined 
to the middle of the wire. Modes that have nodes in the 
middle of the wire, i.e. the second, fourth etc. mode, see 
very little of the potential and go easily through. This 
is especially clear in the case of the longitudinal poten- 
tial were these modes have almost perfect transmission 
as soon as they become propagating. The transverse bar- 
rier is nearly independent of the transverse direction and 
affects all the modes equally. 

The lower panel of Fig ITTl shows the conductance of at- 
tractive potentials. Compared to the symmetric case, the 
elongated longitudinal and transverse attractive barriers 
have higher and lower conductance respectively. The dip 
structure is quite different. The shape of the dips in the 
transverse case are more of the asymmetric Fano form 
than the more common Breit-Wigner form as in the lon- 
gitudinal casej^i In the longitudinal case, we have two 
dips rather than the usual single dip. Because of symme- 
try blocking, there are no dips in the first band as before. 
Even though no dips are seen in the third subband we 
behve they are there. They are just to narrow to be seen 
in the energy resolution of the calculations. 
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Figure 12: (Color online). Conductance of a hard wall quan- 
tum wire in units of Go = 2e^ /h, as a function of energy. The 
scattering potential has the same parameters as the attractive 
potential in Fig. but the center is moved from the middle 
of the wire Xi — L/2 to Xi = 5L/8 and Xi — 3L/4. The width 
of the wire is L = 6 ao and the total number of modes is 
N = 8 SiS before. The inset shows the Vbe~^ contour of the 
scattering potentials. The potentials in the center of the wire 
and the one with Xi — 5L/8 have been shifted by 1 Go and 2 
Go respectively, for clarity. 

The dips are most probably too narrow to be seen in 
experiments. By shifting the center of the potential from 
the middle of the wire, dips in the first subband are no 
longer blocked (cf. Fig. I12II . When the center of the lon- 
gitudinal potential is in Xi = 3L/4 we have two dips that 
are quite pronounced. We will now focus on these two 



dips. 

If we include only one mode, i.e. the propagating mode, 
in the calculations these dips disappear. Both dips reap- 
pear by adding just the first evanescent mode. The higher 
evanescent modes change the conductance curve only 
quantitatively, shifting the dips and affecting their width. 
We therefore deduce that both dips can be explained by a 
coupled mode model restricted to two modes. This kind 
of model as already been put forth to explain the single 
dip in the presence of a single impurityiSfiiSl Nockel and 
Stone noted that a quasi-bound state <I>, satisfying 

2rn dz^ 



+ V22{z) $(z) - [Eq - e2)$(2) (27) 



appears from the evanescent mode. They note that there 
is no need for the dip to be close to the subband edge^l. 
Even though they do not mention it, it is apparent from 
their formaHsm that there is also nothing that prevents 
more than one quasi-bound state to be present, as long 
as there are more than one solution to Eq. Ij27|l with 
Eq — 62 < 0. We therefore suggest that the two dips 
are due to two quasi-bound states originating from the 
first evanescent mode. To confirm this suggestion we 
have solved Eq. Il27p numerically by putting a large box 
around the potential and expanding the bound state $ 
in the plane-wave eigenfunctions of the particle in a box 
with periodic boundary conditions. The energy eigen- 
values obtained are plotted in Fig. showing clearly 
that there are two bound states, that develop into quasi- 
bound states due to coupling to the propagating state. 



■a 



I 

o 




Figure 13: The energy levels of the effective potential V22{z) 
of the longitudinal attractive potential of Fig. 1121 with Xi = 
3I//4. The levels are obtained in a box of length Lbox = 150 



tSQ. 



Depending on the shape and strength of the scattering 
potential the above argument suggests that one could 
see even more dips in each subband. Furthermore, if the 
binding energy fulfills -Eo — £2 < — (^2 — £1), a true bound 
state originating from the second mode would exist. 
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Figure 14: (Color online). The probability density of the scat- 
tering states tpnE hard wall quantum wire in the pres- 
ence of the elongated longitudinal Gaussian potential shown 
in the inset of Fig.lT^shifted to Xi — 3I//4. The incident wave 
is in mode n = 1 and has total energy E = 0.865 (1.084) Ryd 
in the upper (lower) panel. These energies correspond to the 
two dips in the conductance of Fig. 11211 . 



IV. SUMMARY AND DISCUSSION 

In this work we have studied the coherent quantum 
transport in the presence of various reaUstic Gaussian- 
type scattering potentials using the Lippmann-Schwinger 
T-matrix approach. We have calculated the conductance 
and the spatial distribution of the electron probability 
and their dependence on the shape of the scattering po- 
tentials. The formaHsm is quite general and can be ap- 
plied to quantum wires modeled with different confine- 
ment potentials, and in the presence of a general ex- 
tended scattering potential. We have used both hard wall 
confinement and parabolically confined quantum wires. 
The scattering potential has been modeled by Gaussian 
potential, using single and double, symmetric and asym- 
metric Gaussians. By a combination of these we have 
modeled a single impurity, a quantum dot embedded in 
the wire and studied the shape effects of the scatter- 
ing potential. Repulsive impurity smooths out the exact 
quantization of the conductance while attractive impu- 
rity shows the well known dips in the conductance. A res- 
onant transmission appears in the conductance through 
the quantum dot as the energy of the incident electron 
coincides with an energy level in the dot. A longitu- 
dinal barrier affects modes of different symmetry differ- 
ently while a transverse barrier affects all the modes in 
the same way. In an attractive longitudinal potential 
two dips can appear in the same subband. These two 
dips appear due to coupHng of a propagating mode to 
two quasi-bound states that both originate in the next 
evanescent mode. 

Our calculations demonstrate the power and flexibil- 
ity of the current approach. Being able to calculate the 
conductance of a quantum wire in the presence of a gen- 
eral extended scattering potential allows for a multitude 
of systems to be studied. Our Gaussian potentials are 
rather simple but by appropriately choosing the scatter- 
ing potential many different experimental setups can be 
modeled. 



To conflrm our suggested mechanism, in Fig. we 
examine the probability density of the scattering states 
at the energies corresponding to the two dips. The sym- 
metries of the quasi-bound states are as expected. The 
one lower in energy is nodeless in the z direction while 
the other one has a single node. Both have the symmetry 
of the transverse mode n — 2. 

We note in passing that the width of the scattering 
potential can be estimated as the width of the Voe~^ 
contour. The elongated potential in Fig^jhas the width 
0.33 oo in the transverse direction and 3.02 gq in the prop- 
agating z direction. We observe that the quasi-bound 
states are much more extended than one would expect of 
a true bound state in the potential. They also have prob- 
ability density on both sides of the wire. The above indi- 
cates that we are looking at the two quasi-bound states 
from mode n = 2. 



Appendix A: SOLVING THE T-MATRIX 
INTEGRAL EQUATION 

When solving the LS equation for the T-matrix, 
Eq. ltT8|l . numerically, one has to treat the singularities 
of the Green's function with care. By using the fact that 



iS{q~ki) + Siq + ki)) 



kf - g2 



ir; k^ — 2ki 



(Al) 

the integral on the right hand side can be transformed 
into two algebraic terms plus a principal value integral 
{V denotes a principal value) 



dq 



Vnii{k,q)Tin{q,k„ 
kf - g2 



(A2) 



10 



Since 



which further yields 



V 



dq 



kf - 



{A3) 



it is possible to remove the singularity of the integrand 
by a subtraction of a zeroS2i2i. For a general function / 
we have (assuming A; > 0) 



r 



dq 



— oo 
oo 



dq 



/(g) 

kf - q"^ 



'0 

fi-ki) 



dq 



q- 



dq 



f{ki) 



(A4) 



kf 



The value of the integrand at the singularity can, by use 
of I'Hopital's rule, be seen to be ^f'{±k)/2k. Assum- 
ing that /' is nonsingular at ±k, the principal value can 
therefore be skipped as already done. Using this prescrip- 
tion and numerical Gaussian integration for the remain- 
ing integrals, the integral LS equation for the T-matrix 
is transformed into a system of linear equations that can 
be solved by standard linear algebra methodsi^fi. 

The Green's function for the evanescent modes is non 
singular since kf < 0. The integral in Eq. ltT8|l for evanes- 
cent modes is therefor done with numerical Gaussian in- 
tegration without the above special treatment for singu- 
larities. 



Appendix B: SCATTERING STATES IN TERMS 
OF THE T-MATRIX 



/ / dx'dz' — r= 



(B3) 



.X:„(x')V^(r')*+^(r'), 



where we have used the expansion in Eq. JSJ. By the 
definition of the T-matrix, Eq. Ijl6|l , we can thus rewrite 
Eq. as 



nO 
'mE 



(^) + ^ / ' 



\p\e 



~ P 

(B4) 

We have thus expressed the wave function in terms of the 
elements of the T-matrix. By Eq. ljB4|l and the expansion 
in Eq. © we are able to determine the scattering states 
for all r. In the case of propagating modes, the 
singularities of the mode Green's function are treated 
the same way as in App.El For evanescent modes there 
are no singularities and we furthermore have that the 
incoming wave, (f'mEi^)' zero. 



By a spectral representation of the mode Green's func- 
tion, Eq. lO, 



1 /■+°° pip(^-z') 

GIe{^.z')^- dp 



271" ./-oo P^ 



the wave function in Eq. Ijl()|l can be rewritten as 

TTft J^oo '^m ~ P 



(Bl) 



X // dx'dz'e'^P^'xU^')V{r'^:^,E{z')Xn.'{x'), 
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